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the experimental value. We also present power counting arguments that establish a hierarchy for 
electromagnetic one- and two-body currents. 
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I. INTRODUCTION 



The field of nuclear structure and reactions have witnessed continuous progress and re- 
newal of interest with the advent of present and future facilities that are able to provide 
high-intensity beams of very unstable, rare isotopes — so-called exotic nuclei — whose phys- 
ical properties are in the process of being uncovered. They populate areas in the nuclear 
chart far from the valley of stability, resembling more weakly-bound clusters rather than 
a tighter, shell-like structure. A subset in this exotic zoo comprises halo nuclei, systems 
formed by a tightly bound core nucleus surrounded by one or more loosely bound nucleons, 
with a slowly decreasing wave function tail that extends much farther than the effective core- 
nucleon interaction. Such extended and dilute configuration leads to threshold phenomena 
with consequences to low-energy nuclear astrophysics [li |2j . 

An example of astrophysically relevant halo nucleus is boron-8, with dominant configu- 
ration of a beryllium-7 core loosely holding a proton by about a tenth of MeV. This nucleus 
plays a decisive role in our current understanding of neutrino physics. Underground detec- 
tors like SNO and Super-K are mainly sensitive to neutrinos released from the /3-decay of 
in the Sun. The capture rate of a proton by beryllium-7 that produces boron-8, ^Be(p, 7)^B, 
is crucial for determining the initial electron neutrino flux that eventually transmutes to 
other neutrino flavors on its way to detection. The reaction rate for conditions in the solar 
core sharply peaks around 20 keV, the Gamow peak, while experimental data remain above 
100 keV. Measurements at lower energies are extremely difficult, mainly due to the Coulomb 
repulsion. The solar neutrino flux is therefore dependent on theoretical extrapolations of 
current data to lower energies. Available theoretical approaches pHTT] that use p-^Be data 
alone, even with scattering information, are not well-constraining. The mirror system n-''Li 
then becomes an important ingredient for benchmarking purposes — the Coulomb force be- 
tween the initial particles is absent and very low (sub-keV) energies 
are available. It also strongly constrains the ^Be(p, 7)^B if mirror symmetry, which ought 
to have its origins in the (accidental) isospin symmetry of QCD, is invoked. Besides, n-^Li 
has its own astrophysical interest, bridging the path to the formation of heavier elements in 
inhomogeneous big-bang nucleosynthesis models [12]. 

The ^Li(n,7)^Li cross section was calculated recently in effective field theory (EFT) for 
halo nuclei at low energy [13]. In this framework, the tight ^Li core, inert and structureless 
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at leading order (LO), the loosely bound neutron, and the external soft photons are the 
relevant degrees of freedom. The main assumption of the approach is a single-particle 
approximation somewhat similar to models like Refs. [HHUl 115) where the valence nucleon 
interacts with the core via a Woods-Saxon potential. In halo EFT, however, a systematic 
and model-independent expansion of observables is achieved through the use of an expansion 
parameter — formed by the ratio of a soft scale Q, associated with the shallowness of the 
valence neutron, and a hard scale A, related to the tightness of the core. Moreover, the 
formalism guarantees unambiguous inclusion of electromagnetic interactions that preserve 
the required symmetry constraints, such as gauge invariance [121 US]. In Ref. [13] only the 
leading El transition from initial s-wave continuum to the p-wave ^Li ground state was 
considered. The authors found the poorly known p-wave effective range to be the leading 
source of uncertainty. In this work we extend the previous one in several directions, namely, 
we include explicitly the capture to the first excited state of ^Li, the contribution from the 
low-lying p-wave resonance at ~ 0.22 MeV, and elaborate on the power counting for the 
two-body currents. With these extra ingredients, one sets the formalism that paves the way 
to handle the more weakly bound p-^Be mirror system. 

The paper is organized as follows. In Sec. |Tl] we develop the basic theory for the inter- 
actions necessary to calculate the capture reaction ^Li(n, 7)^Li. We present the Lagrangian 
for elastic scattering in the s- and p-wave in the n-|-'^Li system and derive the one-body 
(magnetic moment) and two-body currents. We describe how the couplings in the EFT 
Lagrangian are constrained from available data on low-lying bound and resonance states. 



Results of our calculations are shown in Sec. Ill and compared with some available data and 
potential model calculations. We also present a brief discussion on higher-order terms and 
other degrees of freedom that becomes important at energies slightly above the ones consid- 



ered here. Finally, we present our conclusions in Sec. |IV[ A more technical discussion about 
the power counting of two-body currents in the present approach is given in Appendix |A| 



II. FORMALISM 

The low energy theory for ^Li(n, 7)®Li is constructed out of the spin-parity neutron 
and I lithium-7 core. The final ^Li nucleus contains the 2+ and l"*" ground and excited 
states in the bound spectrum. Very close to the n-''Li threshold only the initial s-wave 
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states are relevant. As the energy increases around 200 keV a pronounced p-wave resonance 
in the initial state, identified as a S"^ state, contributes. It is useful then to first list the 
possible initial and final channels for the reaction. Concentrating on just the s- and p- 
waves, we have in the spectroscopic notation "^^^^Lj: the initial s-wave states ^Si and ^5*2, 
the final p-wave channels ^P2, for the ground state and ^Pi, ^Pi channels for the excited 
state, and initial p-wave resonant state ^P^- The 2"*^ ground state of lithium-8 has the 
quantum numbers of both ^P2 and ^P2 states. It is, however, identified with the symmetric 
combination |2+) = (|^-P2) + |^-P2))/"\/2 [16j. The 1^ excited state is primarily dominated by 
the antisymmetric combination |1+) = (|^-Pi) — \^Pi))/\/2 [7]. The 3+ resonance can only 
belong to the ^P^ channel in the present n +^ Li approach. The leading contribution to the 
capture reaction, which comes from the initial ^Si and ^5*2 states to the 2+ ground state, 
proceeds through El transition due to electromagnetic selection rules. There is a small 
contribution from El capture to the excited 1^ state with a branching ratio of 0.106 [T7] . 
Finally, around 200 keV there is the Ml contribution from the S"*" resonance to the 2+ ground 
state. 



In Ref. [13], the authors used the Clebsch-Gordan coefficient matrices Eq. (B6) to project 
the neutron N{x), lithium-7 C{x) fields onto the spin S = 1 and spin S = 2 channels as 
N'^FiC and N^QijC, respectively. We use the same formalism here to construct all the 
relevant initial and final states. 



A. Interaction 

The operators required for the calculation of the capture reaction ''Li(n, 7)^Li can be 
classified into three categories: (a) s-wave initial state elastic scattering Cg, (6) p-wave 
elastic scattering to describe the 3^ resonance, and analytically continued to describe the 
2+ ground and 1"*" excited states Cp, and (c) two-body currents that are not related to the 
elastic channels £em = O'^^ + O^. 

The s-wave interaction Lagrangian is written as [13] 

= g^'\N^F,CyiN^F,C) + g^'\N^ Q,,C)\N^ Q^fi) + ■ ■ ■ , (1) 

where the "• • • " represent higher derivative terms that are suppressed at low energy. At LO 
there is a single coupling (^^^^ in ^5*2 ('^«S'i) spin channel that is fixed from the known s-wave 
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scattering length a§'^ = -3.63 ± 0.05 fm (a^^^ = 0.87 ± 0.07 fm) [E]. 

As discussed in the next subsection, the description of a low-energy p-wave bound/excited 
state or resonance requires two operators at LO 1^ . It is convenient to work in the 
dimer formalism, where the four-fermion contact interaction is replaced by the exchange of 
an auxiliary dimer field The p-wave interaction Lagrangian then reads 



r{ri) _ aM'' 



^fl + Zi^") C;iV^p(^)C + H. c. , (2) 



where M = M„ + M^, neutron mass M„ = 939.6 MeV, "^Li mass Mc = 6535.4 MeV and H.c. 
stands for Hermitian conjugate. P^^'^ are the projectors for the relevant p-wave channels: 
^Pi, '^P2, ^Pi, ^P2, ^Ps indicated by the index t] and explicitly given in Appendix [b] The 
subscript [j] is a single, double or triple tensor indices as appropriate for J = 1, J = 2 and 
J = 3 state, respectively. The two EFT couplings A^^\ h^^^ are proportional to the two 
required operators at LO and are determined from elastic scattering data as we discuss later. 
For the capture through the 3"'" resonance state, we need the magnetic moment couplings 

O*^ = gnf^NN^ (I ■ AT + g,fir,C^ (J . B) C, (3) 

where cr are the Pauli matrices, J are the angular momentum matrices for spin-3/2 particle, 
-B = V X A the magnetic field, fi^ the nuclear magneton, and gn = 2/t„ and gc = 2kc/3 
are the neutron and the ''Li gyromagnetic ratios, respectively. We take k„ = —1.91304 
and Kc = 3.256427 as the corresponding magnetic moments |2T]. In addition, there are 
contributions from two-body currents 

that contribute to Ml transition from initial ^Ps state to final ^P2 and ^P2 states, respectively. 
The tensor Gijkimq is defined in Appendix |Bj The two-body currents are allowed by symmetry 
so they contribute to the Ml capture. Further, the two-body couplings L^^\ L^'^^ also regulate 
the divergences that appear in certain loop diagrams. The power counting for the different 
interactions will be presented in details in the following and in Appendix |X] 



B. EFT couplings 



The couplings in Eqs. ([T]), ^ can be related to elastic scattering data in the s- and 
p-waves when available. Therefore, it is appropriate to match the field theory to the low- 
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energy amplitude written in terms of the effective range expansion (ERE). In principle, one 
could incorporate relativistic corrections in the EFT amplitude and go beyond the ERE |22] ■ 
However, this is not required at the low energies that we are interested in here. The two-body 
current couplings L^^^ are not related to the elastic scattering data, and are thus determined 



from the capture data (see Sec. III). 

The ERE elastic scattering amplitude in the i-th partial wave is written as 

/i cot 5e - ^ + l^^p2 + l^-^pi + ^p2£+l ' ^ J 

with reduced mass /x, and ae, ri, ti, etc., ERE parameters. Each term is assigned a mo- 
mentum scaling, in general, given by dimensional analysis. Since the ERE parameters are 
associated with the short range (high momentum A) nuclear interaction, naively one would 
expect ai ~ A^^+^, ~ A^^^^, t£ ~ A^^^^, and so on. At arbitrarily low momentum p ^ Q 
one can expand the amplitude Ae in a Taylor series around Q/A = 0, a situation where 
the interaction is weak and perturbative. The more interesting situation arises when there 
are shallow bound or virtual states that one wishes to incorporate in the formalism. That 
implies a rearrangement of the perturbative series, which is only possible if at least one 
of the ERE parameters (usually the scattering length) has a different scaling than the one 
assumed by naive dimensional analysis. The s- and p-wave scattering amplitudes at LO 
then read [ig 1201 [221 123] 

. . / X 27r i . . 

lAoip) = 77 r- , 6 

iAi{p) = 1 ^ . (7) 

/i — 1/ai + T^rip'^ — ip-^ 

For a shallow bound or virtual state in s-wave we set 1/ao ~ Q-, then only one single 
operator is needed at LO. For a shallow p-wave the situation is more subtle. First, only 
one particular fine-tuning of the scattering "length", 1/ai ~ <5^A, is enough to produce 
the shallow state [20]. However, not only one but two operators emerge at LO, since the 
effective "range" term ~ Q^A now scales equally as 1/ai at momenta p ~ Q. Second, 

the unitarity term ip^ ~ is in principle of higher order. The p-wave amplitude Ai 
is then suppressed by Q/A relative to ^o- However, for energies close to the resonant 
state there is a cancellation of the leading terms (kinematical fine-tuning [201 El]) that 
makes — 1/ai + rp^ /2 ~ and promotes the unitarity term ip^ to LO. In that region, the 
amplitudes ^ and ([T]) contribute at the same order. 



The elastic scattering amplitude in EFT is calculated from the interactions in Eqs. ([T]), 
(|2| as shown in Fig. [l] We get for the s-wave amplitude 

with renormalization scale A. Comparing Eqs. ([s]) and ([6]) we find [221 123] 

^^^''HA) = (9) 

where rj = 1,2 correspond to spin channels ^Si, ^82, respectively. 
For the p-wave amplitude we have from Fig. [T] 

[hM]2 2 V 2 ['^''''P / ^ ^ 

using the full dimer propagator 



^D'^-\po,p) =— -, (11) 

1 / 3 , vr 



where C = y— 2/ipo + l^p^ /M — ^0+. Matching the EFT amplitude to the p-wave ERE 



expansion determines the coupling pair (A(''), /iW). As in the s-wave, comparing Eqs. ( |10[ ) 
and ([7]) yields 

_2.,A(^)_.^3_,/,(.)^ 



[/l('')]2 2' 

2^ [h(r<)f -2^ ■ ^^^^ 

The index i] identifies the relevant p-wave channels outlined earlier. The EFT couplings A^''^ 
are therefore fixed in terms of the ERE parameters a|^^\ r^'^ and the renormalization 
scale A. 

In the resonance channel we call 77 = 3 and directly use the relations in Eq. (12). The 



parameters a\ and r\ are related to the known resonance position and width, as shown 
in Sec. HTCl 

For the bound channels 1+, 2+, we follow the procedure used in Ref. [13]. It is more 
convenient to work directly with the location of the pole in the dimer propagator at the 



binding momentum 7'^''^ and its residue Z^'^\ The latter is the wave function renormahzation 



constant calculated using Eqs. (11), (12) as 



where E^'') = [7^'')] V(2/i) is the binding energy. The and channels in 2+ share a 
common binding momentum 7^^^) ^ 57.8 MeV. Moreover, the capture cross section is not 
independently sensitive to the effective range parameter ri in these two spin channels [T3] . 
For this reason we use a common effective range parameter \ This coincident ally gives 
the observed ratio 0.82 of spin channel S* = 2 relative to the total El capture to the ground 
state at threshold [25]. We make a similar simplifying assumption for the 1"^ state and use 
a common effective range parameter ^ for both spin channels ^Pi, ^Pi. In the final cross 



section only the combinations 



j^(2+)]2^(2+) ^ 27r_ 



37(2+) 



,(2+)' 



3^(1+) +^(1+)' ^ 

contribute to the 2"*" and 1^ states, respectively. In Ref. [13], •* ~ —1.55 fm~^ from a 
fit to low energy data from Ref. [2S]- 7^^^-* ~ 41.6 MeV from the known 1"*" excited state 
energy. We will determine ^ and ^ in this work from the known El thermal capture 
rates to the 2+ and l"*" states [17J, respectively. 




+ • 



FIG. 1. Initial s-wave A}^""' and p-wave A\ elastic scattering amplitudes. Double line is the '''Li 
propagator, single line the neutron propagator, dashed line the bare dimer propagator. 
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C. Resonance parameters 



The p-wave EFT couplings for the 3^ i^Ps) state can be related to the ERE scattering 

/Q\ /o\ I I 

parameters a\ and through Eq ( 12 ). Since these ERE parameters are not known we first 



determine them from the known resonance energy and width. At very low momentum, the 
phase shift 6{p) in this channel vanishes due to the centrifugal barrier. Near the resonance, 
6{p) increases rapidly through 7r/2 from below. Thus, cot 6{p) has to go through zero from 
above, i.e., it has a negative slope. We write these conditions as [2] 

cot 6 

d 



= 0, and (15) 

= -c < 0, 



such that c > 0. Near the resonance 



cot 6{E) ^ cot 6{Er) + {E- E,) cot' 6{E,) = -{E - Er)c, (16) 

, 27r 1 27T -Tr/2 

j[(p) = ~ — 

/i p cot 6 — ip np{E — Er) + iTr/2 ' 

and we recover the Breit-Wigner form, identifying c = 2/rr. The resonance position and 

width can be related to the ERE scattering parameters. In the center-of-mass (cm.) frame 

we define = 2fiEr, then 



I ^(3) 2 ^(3) 

p3 cot 5 = -r + cot 6{pr) = 7^ + ^ = 0, and (17) 

il' a\ 'pI ^Pr 



cot 5{Er) = ii A cot 5{pr) = -2/V, - ^ ''^'^ - 



dE prdpr ' a(3)^4 2p2 ^^r. 

For the resonance momentum and width one gets (see also Refs. [281 129] ) 



2 2p^ 

^^' = ^3M^' = — ^ (18) 



Solving these equations we find 

We choose pr = +\/2fiEr so that {af\r^{^^) are negative to be consistent with the Wigner 
bound [30]. If we take the experimental central values for Ej. = 0.222 MeV and = 0.031 



MeV in the cm. frame we find (|l/afV^^^^) = (46.38, -547.1) MeV, and from Eq Q 
we are able to determine the EFT couplings. 



To test the robustness of the above procedure we generate "synthetic data" for the res- 
onance phase shifts from a known nuclear interaction. We take a single-particle potential 
given by the Woods-Saxon form [TT], 



V{r) 



-Vo 



1 + exp 



IVso d 

2 r dr 



1 + exp 



(20) 



where the 1/2 factor in the second term comes from the expectation value of the single- 
particle spin-orbit operator in the channel. We use natural units with h = 1 = c. 
In a study by Tombrello [9], the central potential with a depth vq = 26.42 MeV, spin- 
orbit Vso = 0, range Rc = 2.95 fm and diffusiveness Oc = 0.52 fm was used to reproduce 
the resonance energy. A more recent work from Huang et al. [llj uses vq = 34.93 MeV, 
spin-orbit potential depth Vgo = 10 MeV, Rc = 2.391 fm, and = 0.65 fm. The two 
sets of parameters produce nearly identical phase shifts, Fig|2} Plugging Huang's synthetic 
data into Eq. (jTsl) generates Er = 0.228 MeV and = 0.115 MeV. From Eq. ^ one 



gets {\l/a\ 



(3) 11/3 J3). 



) ' 1 



(30.69, —154.3) MeV, which are then used as input to the EFT 



curve shown in Fig|2} Alternatively, one can extract directly the ERE parameters from 
the low-energy behavior of the phase shifts: (|l/a[^^|^/^, r^"*) = (31.02,-157.6) MeV and 



(|iM' 



(3) 11/3 J3)y 



5 ' 1 



(30.84, —158.9) MeV in the case of Tombrello and Huang's parameters. 



respectively. We also show the Breit-Wigner curve using the extracted resonance parameters. 
We see that the procedure outlined allows EFT to reproduce well the synthetic data at low 
energy. 

We further note that the potential model parameters chosen give a resonance width that 



is about three times larger than the experimental value. As we discuss in Sec. IH, it does 
impact the shape of the Ml curve in the capture reaction. To fit the resonance width to 
experimental data would require tuning yet another parameter (besides the potential depth) 
in the Woods-Saxon potential. In the EFT language, two operators are required to describe a 
shallow p-wave resonance at LO. The corresponding couplings (A*^^\ h^^^) are directly related 
to two experimental data: Er and F^. 



III. CAPTURE CROSS SECTION 

The ^Li(n, 7)^Li capture from E = to about 0.15 MeV total cm. energy is almost 
entirely given by the El transition from initial s-waves. The dominant contribution goes 
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FIG. 2. ^P^ phase shifts. Blue dot-dashed curve using Huang's [llj Woods-Saxon parameters, red 
dashed curve using Tombrello's [9j parameters, black solid curve EFT fitted to Huang's "synthetic 
data" , black long-dashed curve Breit-Wigner form. 

directly to the 2+ ground state, with a small fraction (about 10%) going to the l"*" excited 
state. In a recent work [13] the former was calculated explicitly, but the later was only 
taken into account implicitly with the use of the experimental branching ratio. Around 
E ^ 200 keV the resonance enhances the cross section through an Ml transition to the 
ground state. In this section we derive expressions for these three low-energy mechanisms 
and compare our results with potential model calculations and few available experimental 
data. 

The ^Li(n, 7)®Li cross section is calculated in the cm. frame with p (k) the core (photon) 
momentum and k ■ p = cos 6. The incoming momentum p as well as the binding momenta 
for the ground and excited states 'j^''^^ are assumed to be of the order of the low-energy scale, 
i.e., p 7^'') ~ Q. The photon at LO has |fc| = fco ~ (p^ + [y^'^^'^)/{2jj,) and the Mandelstam 
variable s ~ (M„ -|- Mc)^ = M^. The cm. differential cross section is 



da 


1 


\k\ 


d(f)d cos 9 




\P\ 



647r2M2 2p^ 



where 7^^^) k, 57.8 MeV and 7*^^^) k, 41.6 MeV for the 2+ and I'*' final states, respectively. 

The El capture to the ground state at LO proceeds through the diagrams in Fig.[3j The 
photons are minimally coupled to the ''Li nucleus by gauging the core derivatives Vc — >■ 
Vc + ieZ(.A with = 3 the core charge. We quote the final result for capture to the 
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id) 




FIG. 3. Capture reactions ''Li(n, 7)^Li. Wavy lines represent photons, 
state s-wave interaction. 



represents initial 



state from Ref. [13] here for completeness: 



M 



El 



X 



M 
|1+X|2 

1 



2[7( 



2+)l2 



sin^ 6 

' p2 + [^(2+)]2 \ p2 _^ [^(2+)]; 



1/(2) 



2 [7(2")]= 



3 [7(2+)]2+p2 



(22) 



with a = e^/(47r). The amplitude lA^gf^^l is obtained from the above by replacing Oq 



(2) 



.(1) 



The El transition to the ground state is written as 



da 



(2+) 
El 



1 



dcos^ 327rM2 



2/ip 



(23) 



taking the ^Li ground state |2+) as the symmetric combination of final states. The total 
cross section a{p) is obtained from a straightforward integration over the angle 6. We have 
used a common effective range parameter ^ in ^^^^^ as explained earlier, see Ref. [13] for 
details. Fitting to the thermal cross section o"*-^^-* = 40.56 mb [17j gives ■* = —1.47 fm^^, 
which is close to the ^ = —1.55 fm~^ value one gets using the data from Ref [26j. In this 
work we will use the value obtained from the thermal capture. 

The El capture cross section to the 1"*" excited state, a^^ \ comprises the same set of 



diagrams in Fig. ^except the final state dimer, which is in the 1"*" state. The sum over the 
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final polarization states gives an overall factor of 3 instead of 5. Therefore we have 



da 



Y 



(1+) 

El 



M 
1 



/i 



p2 + [^(l+)]2 1 p2 + [^(1+)]^ 



-Y+Y* 



(24) 



(2) 



2[yi+)]3^^ 

3 p2 + [^(l+)]2 + 



p2+[^(l+)]2i|_/^e^i)|2+|_y^ 



('Pi) I 2 
El 



rfcos^ 327rM2 



2/ip 



(25) 



(2+) 



^ 7 



(1+) 



in Eqs. 



where we made the replacements [/i*^^^)]^^*^^^) — )■ [h^^'^^'^Z^^'^\ 7 
The 1 1'*') state is considered the anti-symmetric combination of the final states. As in the 
2+ case, the amplitude in the other channel spin IA^ei'^^I derived from Oq^'* — )■ Oq^^ with 
a common effective range parameter \ From the thermal capture rate a^^^^ = 4.80 mb 
to the 1^ state we get ■* ~ —1.93 fm^^. The El capture cross section to the ground and 
excited state is shown in Fig. |4j We also show the potential model results for comparison. 
The leading uncertainty in potential model results seems to be associated with the poorly 
known effective range r\ that we determine from the thermal capture rate. We also notice 
that the data set that we call ImhofB is more consistent with the 1/v behavior suggested 
by the Blackmon [26] and Lynn [17] data than ImhofA. 

Next, we consider the Ml capture cross section. It proceeds through an initial p-wave 
state, and therefore, suppressed at low energies. In contrast, the El capture takes place via 
initial s-wave states and displays the known 1/p enhancement at low momentum. However, 
the presence of the 3"'" resonance enhances the Ml contribution around the resonance energy, 
making it comparable to El. In Fig. |5]we show the diagrams that make the leading con- 
tributions to the Ml capture. The first two involve the neutron and core magnetic moment 
couplings, and contribute to both and final states. In the third one the magnetic 
photon couples to the charged ''Li core "in flight" or, in a more classical picture, to the 
electromagnetic current generated by the orbital motion of the charged ^Li core. It arises 
from minimal photon coupling and contributes only to the final state. The last diagram 
contains two-body currents in the respective and ^P2 channels. Naively, counting only 
factors of Q/A, the contributions from two-body currents seem to be more important than 
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FIG. 4. Black long-dashed and solid curves are the EFT results for the El capture to the excited 
state and the total El capture, respectively. The shaded area shows the estimated 30% EFT errors 
in the latter. The results of the potential model code CDXS+ [31j using parameters from Ref. |10) 
and Ref. [9] are given respectively by the blue dot-dashed and red dashed curves. The experimental 
points are from Refs. [HI ES [32l [33] . 




(c) [d) 



FIG. 5. The photon is coupled through the magnetic moment to the neutron (single solid line) and 
^Li nucleus (double solid line). Single dashed line is the ground state ^Li 2+ dressed dimer, double 
dashed line is the '^Li 3"*" resonance dressed dimer. 



the ones from the magnetic moments. We will come back to this point in the following. 
Evaluation of these diagrams in the final channel yields 



X Ul{-p)Q,mUc{p) 



27r/i 



3/3 1 



D(3)(//2/i,0) VS 
I{p;X)- 



/l(3)/l{2+) 



with 



^(p;^)- 3 [^(2+)]2+p2 
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(26) 



(27) 



The various tensors above are defined in Appendix Bl For the P2 channel one gets 



iM 



Ml 



X UU-p)QimUcip) 



27r/i 
27rL(i) 



V3 



(28) 



Numerically the gyromagnetic factors are a/3/2(3(7c/2 + (771/2 + 2jj,ZcMn/M^) ~ 1.8 for 
and a/3/2(3(7c/2 — 3(7„/2) ~ 11 for ^P2- The former is of natural size for a dimensionless 
constant — in Eq. (26), the two-body current dominates for a natural L*^^) ~ 1 and the loop 
contribution is subleading. However, the numerical factor is large and enhances the 



loop contribution beyond the estimates of the power counting. Thus in Eq. (28) the loop 
contribution is as important as the two-body current and both enter at leading order. For 
convenience, we keep the loop contribution and two-body magnetic coupling at the same 
order in both the spin channels. The dependence on the renormalization scale A coming 



from the loop function I{p, A), Eqs. (26), (28), is then cancelled by the two-body couplings 
L^^^ and L^'^\ respectively. We write 



^(2) 



27r 



3 1 2nZ^Mn 

^9c + ^Qn + ^2 



A + /3 



(2) 



3 3 



A + /3(i 



(29) 



2" 2^ 

where /3*^*'''s are renormalized two-body parameters with dimensions of mass. The total cross 
section can be written as 



a 



(2+) 
Ml 



1 7a/i (2+)i2|7{2+)i(P^±]7^ 



2\3 



14 3 M2 



X ■ 



+ 



3 1 

^9c + ^9n + 



M2 



ip 



[^(2+)] 



+ /3 



(2) 



3 3 

9c - T,9', 



ip 



+ /3 



(1) 



(30) 



2""^ 2^V [7(2+)]2+p2 

where the proton mass Mp = 938.3 MeV is used. We summed over the final state dimer and 
photon polarizations, and averaged over the initial spin states. The magnetic moment and 
orbital momentum weights are easy to understand if one compares with the non-relativistic 
quantum operator for the Ml transition, which is proportional to 



M2 



L + QcSc + 9nS^ 



N 



(31) 
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and its expectation value between the initial ^P^ state and the final and states, 
respectively. 

All the elastic scattering parameters have been determined. The final expression in 



Eq. (30) depends on two parameters (3^^\ (5^'^\ that we fit to capture data near the resonance. 
The EFT couplings (A(3),/i(3)) 

were matched to the position and width of the resonance in 
the elastic channel. The values of /3*^*''s primarily affect the height of the cross section 
near the resonance, but not its position or width. If one follows the power counting naively 
then only the two-body currents contribute at LO, and the two (3^^^s are correlated. In 
the resummed amplitude, we find a similar behavior in our fits. Thus we use a common 
/3 = = We find /3 = 170 MeV when we fit to data set ImhofB [32] (more consistent 
with the low-energy 1/v behavior observed experimentally) using the experimental 3"^ width. 
Instead, a fit to the same data set but with the 3"'" width extracted from Huang's potential 
model phase shift provides /3 = 83 MeV. Fitting to the data set ImhofA with Huang's 3"'" 
width give /3 = —44 MeV. The results are shown in Fig. [6| [7j Note that the authors assign 
a 20% error to the data sets ImhofA and ImhofB [32]. This means away from the resonance 
where the cross section is small, the errors are also small. This makes the region where the 
resonance contributes the most the least constraining in the fits. 




3=170 MeV 



83 MeV 

p=-44 MeV 

Huang '10 
Tombrello '65 



-■■f ~i\ 

0.4 0.5 



LAB 



(MeV) 



FIG. 6. Ml capture. Black curve with /3 = 170 MeV, resonance energy and width fitted to 
experimental data. Blue dashed curve with /? = 83 MeV, red dot-dashed curve with /3 = —44 
MeV, resonance energy and width fitted to Huang's potential modek Green circles are results from 
Huang et al., blue squares Ml capture using Tombrello's potential model in code CDXS-I-. 

The fit with a wider resonance, Fj. = 0.11 MeV, seems to describe the data better than 
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FIG. 7. Capture cross section including El and Ml transition. We use r\ = —1.47 fm~ , 
^ = —1.93 fm^^. Black solid curve with /3 = 170 MeV and experimental resonance width, 
blue dashed curve (3 = 83 MeV, red dot-dashed curve with (3 = —44 MeV, and Huang's potential 
model resonance width. Green dots are results from Huang et al. The experimental points are 
from Refs. [HI [Ml EH [33] • 

using the experimental width = 0.031 MeV, see Fig. [7j Smaller errorbars near the 
resonance than the estimated 20% ones would be useful to make stronger statements. Higher 
order EFT corrections estimated to be around 30% (see Appendix |A]) could also make the 
fit better even with the narrower experimental width. With these considerations in mind, 
we see that with the current data the LO EFT does not reproduce the width seen in the 
capture. This could be an indication of the limitation of the single-particle approximation to 
describe the capture reaction near the resonance energy as it is towards the higher momenta 
region of the domain of applicability of the low energy EFT. Nonetheless, the same feature is 
observed in the microscopic calculation of Ref. [6]. The potential models get a wider width 
by coincidence and it is not in agreement with the experimental value F^ = 0.031 MeV. 

To expand the range of applicability of the EFT to slightly higher momenta, one needs 
to include other missing degrees of freedom. The | excited state of the ^Li core (''Li*), 
which contributes only to the spin-1 channel, is the first to consider. It can be incorporated 
explicitly in the present halo EFT, since its energy remains close to the core ground state 
(~ 0.5 MeV) and far from the first breakup channel (~ 2.5 MeV). In the n-''Li cm. system, 
the energy required to probe the ''Li* involves momenta of ~ 80 MeV. In the present work, 
the higher incoming momentum considered is ~ 40 MeV, which justifies having the ^Li* 
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"integrated out". An analogous situation is the Delta resonance in chiral perturbation 
theory (the EFT for pions and nucleons) where the Delta can be integrated out of the 
theory when the energy relative to the pion-nucleon threshold is smaller than the nucleon- 
Delta mass splitting [2H [29], [M] . Further, it can also be shown that the virtual contributions 
of ^Li* to the ground state of ^Li in the neutron capture reaction is a subleading effect |35] . 

Contributions from rf-waves to the El transition, which are in principle suppressed by 
a factor of (Q/A)^, may also become relevant with increasing energy p|0], depending on 
the desired accuracy. The next level of sophistication is the inclusion of alpha and triton 
degrees of freedom in a three-cluster treatment, which shall improve the description of the 
Ml capture reaction [71 [8]. 

IV. CONCLUSIONS 

In the present work we extend the previous halo EFT calculation of the ^Li(n, 7)^Li 
capture reaction to include the complete El transition at leading order, as well as the 
leading Ml capture at low energies. We present model-independent results that quantify 
the current uncertainties in nuclear theory in the single-particle approximation, therefore 
serving as a guide to its limitations and also highlighting where more precise experimental 
input is necessary for improvements. 

We include explicitly the El capture from s-waves to the excited state in '^Li that con- 
tributes about 10% to the cross section at very low energy. The new EFT coupling constants 
associated to this process are completely determined by the binding energy of the excited 
state and the El thermal capture rate to the excited state [T7]. For energies below 100 keV, 
our results show the expected 1/v behavior also seen in potential models, however, differing 
by a sizable overall normalization, directly related to the effective range in the ground state 
channel r\ . This is the dominant source of uncertainty in this low energy region — it could 
be determined from accurate elastic scattering data and partial wave analysis, but due to 
the present lacking of the latter it is poorly known. A fit to El thermal capture data gives 
rf''^ ^ -1.47 fm-^ see also Ref. [I3]. 

The Ml capture proceeds via the S"*" resonance near E ~ 0.2 MeV in the ^Ps channel. 
It is suppressed at very low momentum due to the p-wave initial state. However, near 
the resonance there is an enhanced contribution to the total cross section that needs to 
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be considered. In the halo EFT approach, we include and estimate the size of one-body 
(magnetic moment and orbital momentum couplings) and two-body currents that enters in 
the calculation. The one-body current contribution is consistent with the effective magnetic 
dipole operator used in potential models, but the two-body currents are new ingredients. 
These latter also renormalize loop contributions, and the corresponding couplings (3^^'^ and 
are our free parameters to be determined from the capture data. The power counting 
for two-body currents was studied, and found to contribute at LO and N^LO in the Ml and 
El capture, respectively. 

The available capture data near the resonance have large errors that hampers the quality 
of the fit. However, given the current data near the resonance, it seems that using the 
experimental resonance width Tj. = 31 keV in the S"*" elastic amplitude does a poor job in 
describing the Ml capture data. This was also observed in the microscopic calculation of 
Ref. 0. A significantly larger width, about three times the experimental one, provides much 
better fits. This is roughly the width that one gets in potential models that are tuned to 
the resonance energy. This is a coincidence since, in principle, one should tune the potential 
models to reproduce not only the resonance energy but also the resonance width 
accurately. In the EFT formalism, the p-wave resonance requires two operators at LO that 
can be fixed by matching to the resonance energy and width. The fact that halo EFT is able 



to describe the resonance scattering (as shown in Section II C) but not to reproduce the Ml 
capture might indicate the limitation of the current approach. As discussed in Section |III[ 
the Ml capture could be on the outer edge of the range of applicability. To expand this 
EFT range, other degrees of freedom have to be incorporated. Within the present two-body 
treatment, the inclusion of the ^ excited state in the "^Li core is the first step towards this 
goal. The next, more radical extension is to take the leading configuration of the ^Li as a 
bound state of elementary alpha and triton "cores". In such a three-body approach, not 
only the ground and the | excited states in "^Li could be considered, but also the | state 
which, according to microscopic approaches [3, El [H] , is important to properly describe the 
resonance. Nevertheless, the apparent discrepancy in the input related to the S"*" resonance 
width that is used in the ^P^ elastic scattering and Ml capture reaction is unlikely to be 
resolved with current experimental information. More precise capture data around 0.22 
MeV (where the Ml capture dominates) is needed to conclusively state if the single-particle 
approximation is sufficient to describe the Ml capture in ''Li(?7,, 7)^Li. 
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Appendix A: More on Power Counting 

We have seen earlier that the contribution of the two-body currents relative to the loop 
diagrams in Ml capture scales as L^^^ {gii^Qh^^^h^"^^), where g^i] is a combination of the 
gyromagnetic ratios Qc, Qn that is of natural size in the channel but large in the 
channel. To estimate the size of two-body currents we also need estimates for h^'^''^ that is 
related to the effective range. In this section we present the details of the power counting 
starting with the p-wave elastic channel. This also allows one to estimate the expansion 
parameter Q/A. 

In Sec. njwe have fitted the ERE parameters a^i \ r^'^ in their respective p-channels. We 
found |l/a[''Y^^ ~ 30 - 50 MeV and r'^'^ /2 ~ 100 - 250 MeV. This is consistent with the 
situation in Ref. ^Q\, where only l/a[^^ is fine-tuned to scale as AQ^ while rf'^ ~ A obeys 



the naive expectation. Then in the EPT matching conditions, Eq. ( 12 ) 



-r™ + -A, 



2 ' 2' 

the renormalization scale A ~ Q from the loop momentum is a higher order contribution in 
a Q/A expansion. We get [h^^^^ = -2n/r^C^ ~ 1/A and ^A^*?) = -I / {r^^'^ a\''^) ~ . We 
expand the dimer propagator as 

D(po,p) =/^-2(Po,p) + D_i{po,p) + ■ ■ ■ , (Al) 
and the LO term is given by the bare propagator, 
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since loop contributions enter at higher orders. The subscript indicates the scahng with 
powers of Q and we have suppressed the superscript t] here and in the rest of this section. 
To renormahze the loop expansion systematically we write the couplings as 



h =ho + hi-\ , 

A =A2 + A3 + ■ ■ ■ , (A2) 
stressing again that the subscripts bookkeep the powers in Q. Matching to the ERE we get 

A2 = A3 = -^A=^ + ^A. (A3) 

jjTiai IjjLri firfai 

At next-to-leading order (NLO) the dimer propagator gets a contribution from the one-loop 

self-energy diagram, 

D-iipo,p) =[zD_2(po,p)]^S(po,p), 

.E(po,p) =^^/i(po,p) + ^A3 = - ^^Xe + ifx' + 

ZTT Anfi All 



•C .3C'A 3A 

+ (A4) 
ri 2/iri firfai 



which is a A-dependent result. However, the elastic amplitude in Eq. (10) 



A* 

^ - hl^^D^,{py2ii, 0) - hl^^D^,{py2ii, 0) 
/i /i 



D_i(pV2;U,0)/D_2(pV2/^,0) + ^ 

tlQ 



^+ (A5) 



/i -1/ai + |rip2 \^ -1/ai + ^rip^ 
is A-independent as expected for a physical observable. This is in agreement with the 
expansion of the ERE amplitude up to NLO in Q/A. At LO the binding momentum is 
given by 7 = a/— 2/(riai). The NLO term introduces a double pole which, at least formally, 
is suppressed in the Q/A expansion. To treat the bound state consistently we rewrite the 
p-wave ERE expansion |22| [37] as: 



cot 6{p) =7=^ + + + . . . , 

-1/ai =7^ + ^572 + ■ ■ ■ , 

ri =s + ■ ■ ■ , (A6) 
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and define the EFT couplings as 

2/i 2/iri 2/iri /iri 

The amplitude is now written as 

27r 2p2 



1 + 2^^^^^^"^^ 



ri(7 - ip) 



(A8) 



/i ri(7 + ij9)(7 - ip) 

where the NLO correction contributes a factor of — 37/ri to the residue at the pole p = ij 
without introducing any spurious double pole. This correction to the residue at the pole is 
consistent with the wave function renormalization calculated earlier, 

h^Z = —. (A9) 

37 + ri ^ ^ 

We keep the complete result instead of expanding. A convenient approach to recover the 
complete result at NLO without resumation is to use the so called "zed" -parameterization [55] . 
In this approach one would define, for example, = — 27r(2 ~ 1)/ (37), hi = 3/iqA/ (Att) and 
recover, for the wave function renormalization, 

Z = l + (Z-l) + + + ---, (AlO) 

where Z — 1 ~ Q is treated as a perturbation, see Ref. [5B] for details. Here we will simply 
use the resummed result for the 2+ bound state and the 1"*" excited state. 

So far we have discussed the fine-tuning required to reproduce a shallow (p-wave) bound 
or virtual state. The power counting above also applies to low lying resonances. However, for 
these cases there is an additional fine-tuning that is purely kinematical that was discussed 
in Refs. [20l l2l] . This second fine-tuning requires the loop contribution to be treated non- 
perturbatively at energies near the resonance. Since we consider energies near the resonance 
in our capture calculations, we resum the loop contributions in the 3"*" initial state. 

Next we come back to the role of two-body currents L/lg^qQh"^) in the Ml capture. Since 
/i^ ~ 1/A as shown above, the relative contribution for a natural two-body current scales as 
A/{g[i]Q). At LO and from the specific numerical values of g[i\, one notices that in the 
channel only the two-body current enters, while in the case both two-body and one-body 
(magnetic moment) currents contribute. In principle, for a systematic treatment one could 
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write P^"^^ = /^g^"* + P^i'"' + ■ ■ ■ and perform the perturbative renormalization outlined above 
for the channel, while keeping the full loop contribution in the channel. We verified 
that such a treatment satisfies the power counting estimates. As mentioned in Sec. Ill, we 
keep the loop contribution at LO in both and channels for convenience. 

The scaling of two-body currents that appear in the El case is different than in the 
Ml capture. To keep the discussion fairly general, let us introduce a dimer field 7r(*-™™) 
for the two initial state s-wave channels ^5*2, '^S'l. Then the relative contribution of the 
two-body to one-body current in the El capture is generically LeiA;o/^^*-'^~™'^™V[^^*"''~™^^'^'']) 
where we considered the operator eL^i(j)\jExT^yzTxyzij for transition from to ^P2 ground 
state, for illustration. /^C^-wave) ^ jg ^j-^g vr-nucleon-core coupling and /^i^-^'^^^) the dimer 
propagator in the s-wave. In the power counting one either takes A^''^^^^^-' ~ A to treat s- 
wave interaction as perturbative as would be the case for small natural sized scattering length 
(the '^Si channel for momenta p < 227 MeV) or take A^''"^^^*') ~ Q to treat s-wave interaction 
as non-perturbative as would be the case for large unnatural sized scattering length (the ^5*2 
channel around p ~ 54 MeV). For a natural Lei given by dimensional analysis, the relative 
contribution LEifcoAi^*-'^~™'^™V[^^*'^~™^™^] scales as either for perturbative or as Q for 
non-perturbative s-wave interaction in the initial state. The former is a N^LO contribution 
whereas the latter is a NLO contribution. 



Appendix B: Normalization of states 

In this work we adopt the following definitions for the nucleon and core field operators 
$^f/<,"(p)f/i'"(p) = 2Af„^x'''x''" = 2M«/2,2, (Bl) 

s s 
r r 

{^7v,a(a;),V^jv,fe(2/)} = {fe(aj),V^c7,b(2/)} = S'^'^\x - y)Sab , 
{NM.Nliq)} = {CM.Cliq)} = {27t)'6^'\p - q)6a, . 
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X^*^ and ^^''^ are nucleon and core spinors in the fundamental representations of spins 1/2 
and 3/2, respectively. One-neutron states are defined as 

\p, s) = Nl{p)\0) , ^ (p, s\q, s') = 2Mjv {27rf S^^\p - q)S,,, , 



(Tq 



^3 

S' 



Y.U^\q)e^'^-{N,{q),Nl{pm = ^^-U^^{pm . (B3) 



{2^f 

Analogously for one-core states. Generalization to multi-particle states is straightforward. 
1. Projection operators 

For each partial wave we construct the corresponding projection operators from the rel- 
ative core-nucleon velocity, the spin-1/2 Pauli matrices CTj's, and the following spin-1/2 to 
spin- 3/ 2 transition matrices 

s -— ( 1 \ ^ / v^oi \ ^_J_/^oioo^ 

^V6y0-10\/3y' ^ v^y0 10\/3y'^v^y0010y' 

(B4) 

which satisfy 

SiSj = -Sij - -eijkCTk , SjSj = -6ij - jf^^^} + ^^ijkJk^'^^ , (B5) 

where jf^'^^^s are the generators of the spin-3/2. We construct the Clebsch-Gordan coeffi- 



cient matrices 



% \/3 % 

Fi^ ^(T2Si, Qij = --^(J2{(TiSi GjSi), (B6) 



and define the tensors 



1 



Txyzij — 2 i^^xzi^yj ~l~ ^xzj^yi ~l~ ^yzi^xj ~l~ ^yzj^xi^ ) 

~(~ ^jm^ik ~l~ ^ij^km) 
+ ^Imi^iq^jk + ^jq^ik + ^ij^kq) + {^il^jm^kq + ^il^jq^km) 



^ijklmq g j ^ 



^im^jq ~l~ ^kl^iq^jm ) , (B7) 
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that assure the correct number of independent indices for a given total angular momentum. 
The latter have the following properties, 



Rijlm Rjilm Rijml Rlmij i 
Tijklm Tjii^ijji TijjfjYii Ti^Yikij 1 



ijklmq ^jiklmq ^kjilmq ^ikjlmq ^ijkmlq ^ijkqml ^ijklqm ^Imqijk i 
RijxyRxylm Rijlm i RijxyTxyklm TijkxyRxylm Tijklm i 

G abcijkG ijklmn Gabclmn RabxyG xyclmn • (-^^) 



We introduce the photon vector (e^ ), spin-1 (e^), spin-2 {sij) and spin-3 tensor {Sijk) 
polarizations, obeying the following sums [39] , 



E 



E 



pol. pol.avg. pol.avg. 



ijlm 



pol.avg. 



mq 



ijklmq 
7 



(B9) 



All these elements, together with the matrices Fi and Qij from Eq. (B6), allow us to build 
the operators, in coordinate space. 



2 M M„, 



ixyi 

y 



i^r^^ = V^Fj— -— 1 Rxy^,. 



y 



-^jtT^^ - V^Qa^y ( ^ ~ ]^ 1 ^ijkxyz, (BIO) 
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or in the momentum space, 



~ ^3 



~f5 



p('P^) = p^Q T ■■ = v — ^-T •• = pP.('^^) 

Fijk ^ ~ "^^/^ Q xyG xyzijk — P'^^/^Qxy ~^ Gxyzijk — P -Pijk ^ i 0^^^) 

where we used 

McMn Mc-Mn _ _ + Pc-Pn Pc+Pn 

^^McT^' ' = I^cTl^' P = P -rP r^.(B12) 

The projector operators, defined by 

P^") = (B13) 
where [i] representing collectively the indices in a given channel, satisfy 

J2 Tr[p(^)p(^)t] =1. (B14) 

pol.avg. 

This can be straightforwardly verified from the relations above and 

Tr[FiFj] = 5ij , Tr[Q,,Q;j = Rijim ■ (B15) 
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